AD-A056  384  LOCKHEED  MISSILES  AND  SPACE  CO  INC  PALO  ALTO  CALIF  PA— ETC  F/G  8/7 

DOUBLY  ASYMPTOTIC » BOUNDARY-ELEMENT  ANALYSIS  OF  DYNAMIC  SOIL-ST— ETC (U) 
MAR  78  P G UNDERWOOD*  T L GEERS  DNA001-75-C-0294 


UNCLASSIFIED  LMSC/D624828  DNA-4512T  NL 


1 <*  1 
ADA 
066384 

g 

B 

END 

OATE 

FILMED 

8 78 

DOC 

X 

CO 

<£ 


Q 


ijr#  ? 7/ 


DNA  4512T 


DOUBLY  ASYMPTOTIC,  BOUNDRY-ELEMENT 
ANALYSIS  OF  DYNAMIC  SOIL-STRUCTURE 
INTERACTION 


Lockheed  Palo  Alto  Research  Laboratory 

3251  Hanover  Street 

Palo  Alto,  California  94304 


LEVEL 


31  March  1978 


Q_ 

O 

C_3 


Topical  Report  for  Period  1 June  1975—31  March  1978 
CONTRACT  No.  DNA  001 -75-C-0294 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B344075464  Y99QAXSC06137  H2590D. 

D D c- 

Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  D.  C.  20305 

''  7 ^ 


r 


i 

* 


Destroy  this  report  when  it  is  no  longer 
needed.  Do  not  return  to  sender. 


unclassified 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


REPORT  DOCUMENTATION  PAGE 


GOVT  ACCESSION  NO 


Topical  RepWt*  fa mm 
1 Jun«*SP75 — 31  Mar 

6~  PERFqRM|N..  OR.G  SlPQR 


DOUBLY  ASYMPTOTIC,  BOUNDARY- ELEMENT  ANALYSIS  OF 
DYNAMIC  SO  1 L-STRUCTURE  INTERACTION  . 


LMSC/D624828  I 


P.  G.  Underwood 
T.  L.  /Geers 


. pv  -.<m  \ i,  -AN'.*  AT;  N N/*Mf  AND  ADDR' 

Lockheed  Palo  Alto  Research  Laboratory, 
3251  Hanover  Street 
Palo  Alto,  California 


Sub  td«*lY99QAXSjC06 1-37 


94304 


Director 

Defense  Nuclear  Agency 
Washington,  D.C.  20305 


UNCLASSIFIED 


if  tin  ■*  Acf>uJL 


Approved  for  public  release;  distribution  unlimited 


This  work  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS  Code 
B344075464  Y99QAXSC06 137  H2590D. 


Soil-Structure  Interaction 
Doubly  Asymptotic  Approximations 
Roundnry  Element  Techniques 


This  report  describes  a doubly  asymptotic  (DA),  boundary-element  (BE)  treat- 
ment of  a surrounding  soil  medium  that  offers  considerable  promise  for 
dynamic  soil -structure  interaction  analysis.  The  soil-structure  interaction 
is  reduced  to  a surface  relationship  that  is  asymptotically  exact  at  both 
high  and  low  frequencies.  Governing  equations  for  linear  problems  are 
developed  in  matrix  form  for  application  to  complex  structures.  Numerical 
results  are  presented  for  a two-dimensional  problem  for  which  analytical  
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SUMMARY 


This  report  describes  a doubly-asymptotic  (DA),  boundary-element  (BE)  treatment  of 
a surrounding  soil  medium  that  offers  considerable  promise  for  dynamic  soil— structure 
interaction  analysis.  The  soil— structure  interaction  is  reduced  to  a surface  relation- 
ship that  is  asymptotically  exact  at  both  high  and  low  frequencies.  Governing  equations 
linear  problems  are  developed  in  matrix  form  for  application  to  complex  structures. 
Numerical  results  are  presented  for  a two-dimensional  problem  for  which  analytical  solu- 
tions have  appeared  in  the  literature.  Good  agreement  between  the  DA/BE  and  analytical 
results  is  observed. 
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SECTION  I 


INTRODUCTION 

The  treatment  of  soil-structure  interaction  is  of  considerable  importance  in  analy- 
ses of  the  integrity  of  structures  in  ground-shock  environments.  There  are  currently 
three  basic  approaches  to  the  linear  treatment  of  this  problem:  analytical,  lumped- 
element,  and  finite-element.  Analytical  approaches  are  restricted  to  very  simple  geo- 
metries; hence,  the  results  are  useful  for  providing  insight  into  the  physics  of  the 
problem,  but  the  extension  to  complex  geometries  is  difficult.  The  lumped-element  ap- 
proach, in  which  the  soil  characteristics  are  represented  by  discrete  masses,  springs 
and  dashpots,  is  economical,  but  the  representation  of  actual  soil  behavior  is  crude. 

The  finite-element  (FE)  approach  can  model  the  problem  to  almost  any  accuracy  desired, 
but  the  large  number  of  elements  required  precludes  efficient  computation.  An  approach 
to  achieve  a more  versatile  and  more  economical  method  for  the  treatment  of  these  prob- 
lems would  combine  the  best  features  of  the  different  techniques.  Such  an  approach  is 
pursued  in  this  study:  an  analytical  approximation  of  the  soil-structure  interaction 
is  combined  with  the  modeling  capabilities  of  the  FE  method,  wtiile  avoiding  the  burden 
of  many  elements  in  the  soil. 

This  report  examines  a boundary-element  (BE)  treatment  of  the  surrounding  soil  that 
offers  considerable  promise  for  complex  soil-structure  dynamic  analysis.  The  structure 
is  modeled  through  the  use  of  an  available  FE  code,  and  the  soil-structure  interaction 
is  reduced  to  a surface  relationship  through  the  use  of  a doubly  asymptotic  approxima- 
tion (DAA)  [1],  which  requires  the  application  of  BE  techniques  [2].  The  present  study 

i 

focuses  on  the  two-dimensional  plane-strain  response  of  structures  surrounded  by  an  in- 
finite elastic  medium;  the  extension  to  more  general  problems  is  discussed. 

The  report  first  addresses  the  development  of  the  method:  the  matrix  equation  of 
motion  for  a structure  embedded  in  an  elastic  medium  is  given,  the  doubly  asymptotic 
surface  relationship  is  presented,  and  the  response  equation  for  the  embedded  structure 
is  synthesized.  Then  the  solution  procedure  is  discussed,  and  three  numerical  examples 
are  considered  that  illustrate  the  validity  and  accuracy  of  the  approach. 
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SECTION  II 


GOVERNING  EQUATIONS 

In  this  section,  governing  equations  for  a finite-element  (FE)  model  of  a structure 
and  a boundary-element  (BE)  formulation  of  a first-order  doubly  asymptotic  approximation 
(DAA)  for  the  soil-structure  interaction  are  provided.  These  equations  are  then  com- 
bined to  form  the  response  equation  for  the  embedded  structure.  Finally,  computational 
procedures  for  the  solution  of  the  response  equation  are  discussed. 


2.1  STRUCTURAL  MODEL 

The  matrix  FE  equation  of  motion  for  a linear  structure  embedded  in  an  elastic  med- 
ium through  which  an  incident  disturbance  propagates  is 

£!s  a + Kg  3.  = - (fj  + fs>  (1) 

where  M and  K are  the  mass  and  stiffness  matrices  for  the  structure,  q is  the  struc- 
~s  ~s  ’ ^ 

tural  displacement  vector,  _f  and  _f  are  surface-force  vectors  associated  with  the  in- 
cident and  scattered  waves,  respectively,  and  a dot  denotes  temporal  differentiation. 

The  mass  and  stiffness  matrices  are  easily  obtained  from  any  available  FE  code.  The 
applied  load  is  considered  separable  into  an  incident-wave  force  that  would  exist  if  the 
structure  were  absent  (hence  a known  quantity),  and  a scattered-wave  force  due  to  the 
presence  of  the  structure.  The  scattered-wave  force  constitutes  a troublesome  unknown; 
hence  an  approximation  is  introduced  for  its  evaluation. 

2.2  DOUBLY  ASYMPTOTIC  APPROXIMATION 

A first-order  DAA  is  introduced  to  evaluate  the  scattered- wave  force  f^  [1],  This 
approximation  is  a surface  interaction  approximation,  replacing  the  infinite  volume  of 
external  medium  by  a surface  coincident  with  the  external  surface  of  the  structure. 

The  approximation  is  asymptotically  valid  at  both  high  and  low  frequencies,  as  are  the 
previously  developed  approximations  for  fluid-structure  interaction  [3,4,5]. 


The  development  of  a first-order  DAA  for  linear  soil-structure  interaction  proceeds 
as  follows.  At  high  frequencies,  each  surface  element  of  the  discretized  structure  acts 
as  an  infinite  flat  plate  radiating  plane  waves  into  the  medium.  This  can  be  visualized 
by  considering  that,  for  a fixed  surface-vibration  pattern  oscillating  at  high  frequen- 
cies, the  characteristic  propagation  wave  lengths  in  the  medium  are  short  compared  with 
the  characteristic  wavelength  of  the  surface-vibration  pattern.  For  normal  and  tangen- 


element , 

this  model  y: 

n 

. n 

gSi 

P cd  a4  u4 

t 

. t 

8Si 

P Cs  ai  ui 

(2) 
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where  p is  the  mass  densitv  of  the  medium,  a.  is  the  surf  .e  area  of  the  element,  c,  and 

l d 

c are  the  sound  speeds  for  dilatational  and  shear  waves  in  the  medium,  respectively,  and 
• ^ • t 

u"  and  are  the  normal  and  tangential  scattered  velocities  at  the  surface  of  the  ele- 
ment; see,  e.g.,  [6].  For  an  assemblage  of  elements,  (2)  lead  to  the  matrix  relation 


p A C u' 
~ •v'm  — S 


(3) 


where  A is  a diagonal  element-area  matrix,  is  a diagonal  sound-speed  matrix  for  the 
medium,  and  u'  is  the  computational  scattered-velocity  vector  for  the  surface  elements 
expressed  in  local  coordinates.  Upon  assembly,  the  local  coordinates  in  (3)  are  trans- 
formed to  the  global  coordinates  for  the  problem  as 


S 


C^s 


(4) 


From  (3),  it  is  clear  that  the  external  medium  appears  to  the  structure  as  an  added  dam- 
per in  the  high-frequency  limit. 


Low-frequency  behavior  of  the  medium  is  described  by  the  quasi-static  surface  rela- 


t ion 


“S  ~m  -S 


(5) 


in  which  K is  a surface  stiffness  matrix  for  the  medium.  In  this  limit,  die  external 
'm 

medium  appears  to  the  structure  as  an  added  stiffness  embodied  in  K . The  construction 

~m 

of  K is  discussed  in  Section  2.4. 

To  construct  the  first-order  DAA,  (3)  and  (5)  are  added  to  obtain 


v = P G AC  G u + K u 
— S ~ — - m ~ — S ~in  — S 

_T  , 


(6) 


where  the  transformation  of  (3)  as  = G g'  results  from  (4)  and  the  fact  that  virtual 

T T 

work  must  be  independent  of  the  coordinate  system  used,  i.e.,  (5ij)  g = (Ju‘)  g'. 

It  is  easy  to  see  the  doubly  asymptotic  nature  of  the  surface  approximation.  At  low 
frequencies,  the  velocity  vector  is  small  relative  to  the  displacement  vector,  so  that 
the  scattered  force  is  essentially  given  by  the  static  stiffness  relationship;  at  high 
frequencies,  the  reverse  is  true,  so  that  the  scattered  force  is  essentially  given  by 
the  radiation  damping  relationship.  In  the  intermediate  frequency  range,  the  DAA  is,  of 
course,  in  error;  the  purpose  of  the  numerical  results  presented  herein  is  to  indicate 
the  magnitude  of  that  error.  If  numerical  calculations  demonstrate  the  need  for  an  im- 
proved approximation,  one  may  be  derived;  for  fluid-structure  interaction,  an  improved 
DAA  has  been  developed  that  substantially  outperforms  the  original  [5]. 
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2.3  RESPONSE  EQUATION 


In  linear  problems,  not  only  is  the  surface-force  vector  separable  into  incident- 
wave  and  scattered-wave  components  [see  (1)],  but  the  surface  displacement  vector  u is 
also  separable  such  that  11  = ii^  + u^.  Hence  (1)  and  (6)  may  be  combined  to  give  the 
doubly-asymptotic  response  equation  for  the  embedded  structure 

M q + pDTG^ A C G D q + (K  + I)1  K D)q  = f + pD?GTA  f C D u + DTK  ut  (7) 

s ~ ~ ' m~  - -1  s - ~m~  -I  m 1 m— I 

T 

where  u and  f have  been  transformed  as  u = P q and  f = D g , i.e.,  D selects  the  struc- 
tural  degrees  of  freedom  that  define  the  soil-structure  interface.  In  (7),  M and  K 
are  readily  provided  by  an  FE  structural  analysis  code,  DC  A C^G  D is  easily  computed, 
f.j  and  Uj  are  known,  and  is  determined  through  the  application  of  boundary-integral- 
equation  techniques  as  now  described. 

2.4  MEDIUM  STIFFNESS  MATRIX 


The  basic  boundary-integral  equation  in  two  dimensions  is  [2] 


uk(P) 


/ 


Tk' (P,Q) 


u*(Q) 


dL(Q)  = 


I 


Uk?(P,Q) 


t f (0) 


dL(Q) 


(8) 


k k 

where  P and  Q are  surface  points,  u and  t are  surface  displacements  and  tractions, 
k£  ki 

respectively,  T and  U are  Green's  functions  for  the  boundary,  and  k = 1,2  and  i = 
1,2  are  the  Cartesian  coordinate  indices.  Through  division  of  the  structure's  (two- 
dimensional)  external  surface  into  a series  of  boundary  elements,  (8)  may  be  expressed 
in  matrix  notation  as 


S u = Ft 

in  which  the  2x2  elements  of  S and  F are  given  by 


(9) 


„ki 

3 . . 

1J 


= + 
ij  ki 


Fk*  = 
1J 


( Tk!  C!  dL. 
/ i J J J 

Lj 

f UkJ-  C dL 

Ju  3 3 3 


(10) 


where  5 and  5,  are  Kronecker  deltas,  i and  [ are  boundarv-element  indices,  C.  and 
ij  k7  .1 

5.  are  assumed  BE  shape-functions,  and  L.  is  the  length  of  the  jt_h  boundary  element.  For 

J 1 , ki  „kH  . , r ~ 

the  two-dimensional  plane-strain  case,  the  kernels  T..  and  l are  given  by  [z,/j 


C r 3 r 

k*  _ 13_  ii 

ij  r . . Sin 
ij  L .) 


C4+2rl),t  + C4< 


k e " 

n.  r..  „ -n.  r . . ,) 
.1  tj,f.  j ij,k 


ki 

U . = C.  (6  C,  In  r - r.  r , . „) 

ij  1 ki  2 ij  ij,k  i j , P. 


(11) 
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where  C^,  C^,  and  are  material  constants,  r^  is  the  distance  from  a node  point 

on  the  ith  element  to  the  variable  (field)  point  of  integration  on  the  j th  element,  n. 

k J 

is  the  unit  normal  to  the  surface  of  the  j^h  element,  n^  is  the  cosine  of  the  angle 
between  n^  and  the  kth  Cartesian  direction,  and  a subscript  following  a comma  repre- 
sents spatial  differentiation  with  respect  to  the  indicated  Cartesian  coordinate  at 
point  3.  In  the  present  implementation,  the  displacement  and  traction  shape-functions 
and  are  assumed  to  be  constant  over  the  jth  element,  so  they  may  be  brought  out 

from  under  the  integral  signs  in  (10).  The  numerical  techniques  used  to  evaluate  the 
integrals  in  (10)  are  discussed  in  the  appendix. 

Once  the  matrices  in  (9)  have  been  generated,  it  is  a simple  matter  to  obtain  the 
medium  stiffness  matrix;  because  F is  nonsingular,  it  can  be  factored  to  obtain 

t = F_1  S u = K u (12) 

— ~ ~~  — ~m  — 

As  the  preceding  development  is  not  based  on  variational  principles,  the  derived 

stiffness  matrix  may  not  be  symmetric;  therefore  K is  symmetrized  before  it  is  used  in 

~m 

(7). 

The  brevity  of  the  preceding  BE  formulation  is  appropriate,  in  view  of  the  exten- 
sive coverage  of  the  subject  provided  in  [2].  The  emphasis  here  has  been  on  the  spe- 
cific approach  of  this  study;  it  has  been  found  to  be  most  economical,  especially  the 
use  of  numerical  integration  to  evaluate  the  matrix  elements  defined  in  (10).  The  appar- 
ently new  technique  of  using  boundary  integral  equations  to  define  a medium  stiffness 
matrix  is  valuable,  in  that  it  facilitates  the  use  of  the  form  (12).  This  form  is 

required  for  an  efficient  marriage  of  an  FE  structural  model  and  a BE  soil  model.  There 
are  also  improved  forms  of  the  BE  method  available  that  utilize  higher  order  shape-func- 
tions to  describe  boundary  displacements  and  tractions,  as  well  as  sophisticated  isopara- 
metric-element  representations;  these  procedures  are  reviewed  by  Cruse  [8].  The  simple 
approach  used  in  this  study  is,  however,  adequate  for  the  purposes  of  the  present  investi- 
gation. If  software  were  to  be  constructed  for  production  analysis,  the  incorporation  of 
refined  BE  techniques  would  be  appropriate. 

2.5  SOLUTION  PROCEDURE 

The  doubly  asymptotic  response  equation  for  the  embedded  structure,  (7),  has  the 
form  of  the  standard  matrix  equation  of  structural  dynamics;  hence,  the  solution  of  (7) 
may  be  accomplished  with  well-established  techniques.  For  the  linear  response  problems 
considered  here,  the  integration  of  (7)  is  performed  in  accordance  with  the  trapezoidal 
rule  [9).  The  equation  solver  used  with  the  time  integrator  is  the  skyline  format  pro- 
cedure of  Felippa  [10]. 
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A study  of  (7)  on  a terra-by-term  basis  is  informative.  The  mass  matrix  produced 
by  REXBAT  [11],  the  structural  finite-element  code  used  in  this  study,  is  diagonal;  a 
consistent  mass  matrix  could  be  used,  however,  without  unduly  complicating  the  solution. 
The  damping  matrix  is  highly  banded  in  all  cases  and  presents  no  computational  diffi- 
culties. The  stiffness  matrix,  on  the  other  hand,  may  be  nearly  full,  due  to  the  added 
stiffness  terms.  (Note  that  the  matrices  generated  from  (10)  are  full.)  For  the  simple 
examples  considered  here,  this  dense  stiffness  matrix  presents  no  difficulty.  However, 
for  large  systems,  the  compact  bandwidth  (low  connectivity)  of  the  structural  model, 
which  is  needed  for  efficient  solution,  would  be  lost  through  the  addition  of  the  fully 
populated  medium  stiffness  matrix.  To  overcome  this  problem,  a staggered-solution  ap- 
proach, such  as  the  one  developed  for  fluid-structure  interaction  analysis  by  Park, et 
al. [12]  should  be  considered  for  large  systems  of  equations.  The  forcing  function,  i.e. 
the  right  side  of  (7),  may  look  complicated,  but  each  term  is  known  and  the  load  vector 
is  easily  computed  by  simple  matrix-vector  multiplication  and  vector  addition. 

The  power  of  the  present  approach  is  certainly  evident  for  engineering  applications 
as  the  FE  and  BE  methods  enjoy  direct  applicability  to  the  complex  geometries  of  engi- 
neering structures.  Furthermore,  the  doubly  asymptotic  response  equation  for  the  em- 
bedded structure  is  merely  the  second-order  ordinary  differential  equation  of  structural 
dynamics . 
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SECTION  III 


NUMERICAL  RESULTS 

In  this  section,  numerical  results  for  the  transverse  excitation  of  an  infinite, 
circular  cylindrical  cavity  and  for  two  infinite,  circular  cylindrical  shells  by  an 
incident  plane,  dilatational  wave  are  compared  with  corresponding  analytical  solu- 
tions. Problem  geometry  and  notation  are  shown  in  Figure  1;  in  all  cases,  plane  strain 
response  is  assumed.  The  coincident  finite-element  and  boundary-element  grids  for  all 
three  problems  consist  of  40  elements  of  equal  length.  The  finite-element  shell  models 
incorporate  straight  beam  elements  with  elastic  moduli  modified  for  replication  of  plane 
strain  conditions. 


The  results  are  presented  in  nondimensional  form.  Length  is  normalized  to  a,  time 

2 

is  normalized  to  a/c,,  and  stress  is  normalized  to  pc  = A + 2p,  where  A and  u are  the 
d d 

Lame  coefficients  for  the  medium. 


3.1  INCIDENT  WAVE 

A plane  dilatational  step-wave,  characterized  by  a compressive  pressure  and 
moving  in  the  x ^-direc t ion , can  be  described  in  terms  of  a -tond  intensions  1 scalar  poten- 
tial as 

♦ t P0(t  - x,  - l)2  H(x  - xt  - 1)  (13) 

where  t is  nondimens  anal  time,  x^  is  nondimensional  position  along  the  x^  axis,  and  H 
is  the  Heavyside  operator.  For  this  incident  wave,  the  shear  potential  is  zero  [6]. 


The  incident-wave  force  vector  _fj , which  appears  on  the  right  side  of  (7),  is  ob- 
tained as  follows.  First,  the  elements  of  the  incident-wave  computational  stress  vector 
in  global  (xpX9)  coordinates  are  determined  by  application  of  classical  continuum  form- 
ulas to  (13)  (6):  this  yields 


Ii 

k? 

rli 


- (X  + 2u  6 lk) 


P0  H(T  “ Xli  " U 


= 0 


(14) 


where  x^  denotes  the  x^-position  of  the  ith  surface  mode.  Second,  a global  stress  <*c- 
tor  is  constructed  from  these  elements,  which  is  then  transformed,  on  the  basis  of  Mohr's 
circle,  into  local  coordinates  as  £'  = M £ . Finally,  the  force  vector  in  local  coordin- 
ates is  determined  as  f ' = - A o',  which  is  then  transformed  into  global  coordinates  to 

— I ^ I 

yield  [c . f . , (6) J 


ii 


- U1  A M £x 
11 


(15) 


The  incident-wave  displacement  and  velocity  vectors  and  uT,  which  also  appear  on 

the  right  side  of  (7),  are  obtained  from  the  classical  relation  = 3 <f>/3  x . This 

k 

relation  and  (13)  yield  as  the  elements  of  these  vectors 


UIi  ~ 5 lk  V " xli  " L)  H(l  ' xli  ' 


V = 6 ..  P H(T  - x, . - 1) 

Ii  lk  o li 


(16) 


3.2  CIRCULAR  CAVITY 

The  cavity  problem  is  formulated  simply  by  taking  M = ^ = 0,  which  reduces  (7)  to 
a first-order  equation.  A comparison  between  results  obtained  by  the  present  method  and 
analytical  results  presented  in  [13]  is  provided,  for  step-wave  excitation,  in  Figure  2. 
Minor  discrepancies  exist  between  the  DA/BE  and  analytical  response  histories  at  early 
times.  At  late  times,  both  sets  of  response  histories  approach  the  appropriate  asymp- 
totes [1,  14];  these  asymptotes  are  x-4,  x-1  and  x+2  for  0=0°,  90°,  and  180°,  respect- 
ively. 

3.3  CONCRETE  SHELL  IN  SLOW  GRANITE 

The  second  check  problem,  the  response  of  a concrete  shell  to  an  incident  wave  of 

rectangular  pressure-profile,  is  also  taken  from  [13].  The  nondimensional  parameters 

for  this  problem  are  h/a  = 0.01,  P /p  = 0.865,  c /c,  = 0.63,  c /c,  = 1.87,  v = 0.25,  and 

o s d o cl 

v = 0.2;  the  duration  of  the  incident  rectangular  pulse  is  10.  I)A/BE  and  analytical  dis- 
placement histories  for  this  problem  are  compared  in  Figure  3.  In  this  figure,  as  in 
Figure  2,  the  DA/BE  responses  generally  tend  to  lag  behind  their  analytical  counterparts. 
As  discussed  in  [1],  this  tendency  is  the  result  of  excess  radiation  damping  introduced 
by  the  DAA.  Also  of  interest  is  the  DA/BE  prediction  of  shell  response  at  0=0°  before 
x=1.53,  which  is  the  earliest  time  a disturbance  can  reach  that  point  [15];  this  nonphys- 
ical result  illustrates  that  the  DAA  is  not  a wave  propagation  approximation.  Despite  its 
deficiencies,  however,  the  DAA  produces  results  that  nowhere  differ  from  their  analytical 
counterparts  by  more  than  10%  of  the  peak  response  and  also  approach  the  proper  late-time 

asymptote.  The  latter  characteristic  attests  to  the  correctness  of  the  K -calculation. 

~m 

Although  the  DAA  tends  to  overestimate  radiation  damping,  its  inclusion  is  absolute- 
ly necessary  for  an  accurate  treatment  of  abrupt  soil-structure  interaction.  This  is 
indicated  in  Figure  4,  where  displacement  responses  corresponding  to  the  DA/FE  responses 
of  Figure  3 have  been  computed  from  (7)  with  C set  equal  to  0.  As  one  would  expect, 
the  highly  oscillatory  response  thus  calculated  produces  extremely  poor  stress/strain 
rest]  Its. 
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DISPLACEMENT  (q^a) 
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Figure  2.  Displacement  Response  of  a Circular  Cavity 
to  an  Incident  Step  ;.*ave  [solid  symbols 
denote  results  from  (13)]. 
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Figure  3.  Displacement  Response  of  a Concrete  Shell 
in  Slow  Granite  to  an  Incident  Wave  of 
Rectangular  Pressure-profile  [solid  symbols 
denote  results  from  (13)]. 

L 

15 

3.4  CONCRETE  SHELL  IN  GRANITE 


The  final  check  problem,  the  response  of  a concrete  shell  to  an  incident  step-wave, 
is  taken  from  [1],  The  appropriate  nondimensional  parameters  here  are  h/a  = 0.05,  PQ/P  = 
1.0,  cs/c£j  = lA/3,  CQ/Ccj  = l/\/2,  and  v = = 0.25.  Velocity  response  histories  at  0 = 

0°  i.  180° are  shown  in  Figure  5,  corresponding  to  DA/BE,  DA/analytical  and  exact/analyt- 
ical  treatments  of  the  structure-medium  interaction.  It  is  seen  that  the  DA/BE  and  DA/ 
analytical  results  are  in  almost  perfect  agreement,  which  is  most  reassuring.  Premature 
initial  response  at  points  in  the  shadow  region  and  excessive  radiation  damping  charac- 
terize the  DA  results  here  as  they  did  in  Figure  3.  The  associated  error  is  modest,  how- 
ever, with  all  results  coalescing  at  late  times.  Stress  response  histories  in  the  middle 
and  inner  fibers  of  the  shell  at  6 = 90°  are  shown  in  Figure  6.  Here,  some  minor  dis- 
crepancies between  the  DA/BE  and  DA/analytical  results  appear:  near  t = 0,  the  DA/BE 
histories  exhibit  a more  realistic  delay  before  a stress  response  appears;  near  t = 1.5, 
short-term  reversals  in  stress  appear  in  the  DA/BE  histories,  whereas  the  analytical 
histories  are  smooth;  finally,  at  late  times,  the  DA/BE  asymptotic  stress  values  are 
slightly  less  than  their  analytical  counterparts.  Much  larger  discrepancies  exist  be- 
tween the  DA  results  and  the  exact  results,  especially  during  the  period  4 < t < 12. 

Even  here,  however,  the  error  never  exceeds  15%,  which  is  generally  acceptable  for  engi- 
neering analysis. 
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Figure  5.  Velocity  Response  of  a Concrete  Shell 
in  Granite  to  an  Incident  Step-wave 
(Exact/Analytical , DA/Analytical , DA/BE). 
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STRESS  AT  90°  (o0/dc^) 
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Figure  6.  Stress  Response  in  the  Middle  and  Inner  Fibers  of  a 
Concrete  Shell  in  Granite  to  an  Incident  Step-Wave 
(Exact/.,  ilytical,  DA/Analvtical , and  DA/BE). 
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SECTION  IV 


CONCLUSION 

The  numerical  results  of  the  previous  section  indicate  that  the  doubly  asymptotic 
approximation  of  [1]  offers  considerable  promise  for  the  satisfactory  treatment  of  dyna- 
mic soil-structure  interaction.  In  addition,  the  boundary-element  methods  described  in 
Section  II  and  the  Appendix  constitute  a firm  technology  for  application  to  engineering 
structures  with  complex  surface  geometries. 

The  extension  of  DA/BE  methods  for  the  treatment  of  nonlinear  soil  response  is  dis- 
cussed in  [16].  As  one  would  expect,  this  is  a major  effort,  requiring  consideration  of 
the  state  of  the  medium  at  points  removed  from  the  surface  of  the  structure.  In  spite 
of  the  difficulties,  however,  the  nonlinear  problem  appears  to  be  yielding  to  the  new 
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APPENDIX 

This  appendix  discusses  the  numerical  approach  used  to  evaluate  the  integrals  in 

k£  k£ 

(10)  for  determination  of  the  matrix  elements  S..  and  F 

iJ  ij 

First,  the  boundary  is  divided  into  2-D  boundary  elements,  each  with  a centrally 

k£  k£ 

located  node.  For  a single  calculation  of  S„  and  F..,  fixed  values  are  assigned  to 
i,j,k  and  l,  and  a circle  is  fitted  to  the  nodal  points  j-1,  j and  j+1 ; this  completely 
determines  the  center  and  radius  of  the  arc  describing  the  jth  element.  The  ends  of  the 
j th  element  are  then  point  on  the  arc  half-way  between  points  j-l  and  j and  point 
on  the  arc  half-way  between  points  j and  j+1.  Hence  L.  is  the  arc  length  between 
points  j-%  and  j+^,  and  the  unit  normal  anywhere  on  the  element  is  completely  defined. 

P f 

Second,  the  displacement  and  traction  shape-functions,  and  C.,  are  taken  as  unitv, 

so  the  integrals  in  (10)  only  involve  the  kernals  (11).  In  this  connection,  it  is  import- 

k£  k£ 

ant  to  remember  that  and  lT^  pertain  to  a f ixed  point  (the  nodal  point)  on  the  ith 

element,  but  to  a variable  point  on  the  j_th  element.  For  j ^ i,  the  geometric  quanti- 
ties in  (11)  are  easily  determined  as 

I 2 2 ^ 

r..  = (x,.  - x15)  + (x,;  - x ) 


rij,k 


(x  - x )/r . . 
kj  ki  ij 


i ’ 1 J 

~'T~  = n.  r..  . + n.  r..  0 (17) 

3 J id,2 

and  Simpson’s  rule  is  used  to  evaluate  the  integrals  with  points  j and  as  the 

integration  points.  For  j = i,  special  evaluation  methods  arc-  used,  as  described  in  the 
following  paragraph. 

kJ. 

With  regard  to  the  integral  of  T it  may  be  shown  [7]  that  [see  (11)] 

I r7j  tf.  „ C.  + 2r..  . r,.  .)  dL  = - r(C , + 1)  6, 

#iir»n  k ? 4 1 1 , k 1 1 , P.  i a 'v 


Jr?!  C/  + 2r,.  , r. 

. ii  Dn.  k 9 4 li, k 1 

L . l 

l 

f (n^  r - n!  r.  ) dL. 

Jx  l ix,2  l ll, k l 


where  the  first  i-subscript  of  the  doubly  subscripted  variable  r..  refers  to  the  fixed 
nodal  point  for  the  it_h  element,  and  the  second  i-subscript  of  r..  and  any  single  i - 
subscript  refers  to  a variable  point  on  that  element.  With  regard  to  the  integral  of 


preceding  page  runk-not  filmed 


Hi  ./* 


U^,  it  may  be  shown  that  [see  (11)] 

J hi  dh  ' ri«-«  («"  r1(1-4)  - *)  + '1(i+»  (*»  ii(1+4)  - i)  (19) 

where  the  subscripts  (i— %)  and  (i+%)  refer  to  the  end  points  of  the  ith  element.  Inte- 
nd 

gration  of  the  second  term  in  the  expression  for  U_._^  [see  (11)]  is  performed  by  means  of 

Simpson1  s rule,  with  the  nodal  point  i and  the  end  points  i-%  and  i-f-%  as  integration 

points.  In  this  exercise,  the  second  of  (17)  is  used  directly  to  evaluate  r r 

ii,k  ii,Jc 

at  the  end  points,  while  it  is  used  at  the  nodal  point  in  conjunction  with  a Richardson 
extrapolation  [17]  of  the  form 


ri(i-e),k  ri(i-e),2.  + 2ri(i-2e),k  ri(i-2e),Sl 

+ 2ri (i+2e) ,k  ri(i+2c),i  " ri(i+E),k  ri(i+e),i  (20) 

where  for  example,  denotes  the  value  of  k [obtained  from  the  second  of 

(1^)J  that  pertains  to  the  ith  nodal  point  and  to  a fixed  point  located  between  the  nodal 

points  i-1  and  i at  a distance  e from  nodal  point  i;  here,  e has  been  taken  as  0.05  L . 

i 

k£ 

Finally,  each  value  of  F is  scaled  through  division  by  L . . This  scales  the  trac- 
k . 1J  J 

tions  t so  that  they,  in  effect,  become  nodal  forces,  producing  a stiffness  matrix  K 
J ~m 

of  the  standard  FE  form. 


(n)  (n)  _ 1 

ii,k  ii,J,  2 
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